Low-Lying Magnetic Excitation of the Shastry-Sutherland Model 



K. Totsuka* 

Department of Physics, Kyushu University, 6-10-1 Hakozaki, Fukuoka, 812-8581, Japan 

S. Miyahara, and K. Ueda 
Institute of Solid State Physics, University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba, 277-8581, Japan 

(19 June 2000) 

By using perturbation calculation and numerical diagonalization, low-energy spin dynamics of the 
Shastry-Sutherland model is investigated paying particular attention to the two-particle coherent 
motion. In addition to spin-singlet- and triplet bound states, we find novel branches of coherent 
motion of a bound quintet pair, which are usually unstable because of repulsion. Unusual dispersion 
observed in neutron-scattering measurements are explained by the present theory. The importance 
of the effects of phonon is also pointed out. 
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In recent years, low-dimensional spin systems with 
a spin gap have been a subject of extensive re- 
search. Among them, a two-dimensional antiferromagnet 
SrCu2 (1303)2 is outstanding in its unique features. These 
include, (i) spin-gapped behavior Q|, (ii) magnetization 
plateaus J0| , and (iii) unusual low-energy dynamics || . 
In Ref. Jl|, it was pointed out that SrCu2 (603)2 may 
be modeled by the 5=1/2 Heisenberg model on the 
Shastry-Sutherland (SS) lattice || (Shastry-Sutherland 
model, hereafter). This sparked experimental- and theo- 
retical researches on interesting features of the Shastry- 
Sutherland model @-[uJ]. 

Strong geometrical frustration of the SS model allows 
a simple dimer-product to be the exact ground state [^) . 
In the zeroth-order approximation, a triplet excitation 
above the dimer-singlet ground state is created by pro- 
moting one of the dimer singlets to triplet. In this letter, 
we consider such particle-like excitations and show that 
interesting two-particle motion (bound states) is possible 
because of unusual dynamical properties. 

The unit cell of the SS lattice contains two mutually 
orthogonal dimer bonds (we call them A and B; see Fig.l) 
and the Hamiltonian can be written as a sum of local 
Hamiltonians acting only on either A- or B dimers and 
those acting on both A and B dimers: 

H = 22 H a + 22 7~L a ,f3 ■ (1) 

a=A,B (c«, l 3)=(A,B) 

In terms of the hardcore triplet bosons t on dimer bonds 
p2[ , the local Hamiltonians are given by H a = J[~j + 
t+ -t a ] and n a ,p = %[8ig(a,j3){iil-(ti3 x t a ) + (h.c.)} + 
T Q -Ta], where T a .p denote the S = 1 operators and 
the sign factor sig(a, f3) equals 1 when the arrow on a 
horizontal bond (/3) is emanating from the vertical one 
(a) and —1 otherwise (see Fig.l). 

As is easily seen, a unique geometry of the SS lattice 
allows neither (bare) one-particle hopping (t^ (x)-t(y)) 
nor pair creation/annihilation of triplets; non-trivial one- 



particle (triplet) hopping is generated only perturba- 
tively Q (it occurs at (J'/ J) 6 and higher). 
Correlated hopping 

Although one-particle hopping is strongly suppressed, 
the situation is dramatically different for two-particle 
cases. The 3-point vertices (e.g. t^-(t x t)) contained 
in the Hamiltonian make non-trivial two-particle hop- 
ping like Fig. 2 possible already at (J'/ J) 2 . Note that 
only one of the two particles hops and that the other is at 
rest merely to assist the hopping. Therefore, we call such 
processes correlated hoppings; two triplets close to each 
other can use this new channel of two-particle motion to 
form various bound states. Although the relevance of 
correlated hoppings in the SS model was pointed out in 
H in the context of magnetization process, the effects of 
it would be most highlighted in the low-energy dynamics. 

Usually, correlated motions are only higher-order cor- 
rections to the dominant one-particle processes. How- 
ever, if the one-particle processes are strongly suppressed 
for some reasons, correlated ones would play an impor- 
tant role in the low-energy dynamics. The Shastry- 
Sutherland model is a candidate for such systems. Below 
we demonstrate that the correlated hoppings can explain 
the unusual dispersion observed in inelastic neutron- 
scattering (INS) experiments 0. 
Two-triplet motion 

Collecting all the two-particle processes up to (J'/J) 3 , 
we found that the two-particle coherent motion is closed 
within the four states (Fig. 2) which are decoupled from 
other states where two triplets are far apart. Because of 
this special property, if we take the following four relative 
configurations (Fig. 2) 

|*(R e ,p tf )) = (| (p)),|6(p)),|c(p)),|d(p))) (2) 

as the basis, computation of the two-triplet spectra re- 
duces to diagonalization of a four-by-four (up to 3rd or- 
der) hopping matrix. Although the form of it is the same 
as that given in Ref. Cffl , the elements now are functions 
of spin-1 operators Ti and T2 of the two triplets; for ex- 
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ample, the interactions Vnn and Vnnn between nearest- 
neighbor- and next-nearest-neighbor pair are given by 
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(3) 



The meaning of the hopping amplitudes Jnn and J^a 
can be read off from Fig. 2. Note that p is defined with 
respect to the chemical unit cell differently from that used 
in Refs. [§,1). For completeness, we add another two- 
particle interaction between a 3rd-neighbor pair (| • •• |): 



V 3rc 



{J'f , HJ'f 

2J 4J 2 



Ti-To 



(4) 



which creates immobile (at this order) bound pairs with 
energy 2Vq{J,J') + V 3] -d for S'tot = 0,1. In the pertur- 
bative regime, a pair with total spin S'tot = or 1 feels 
attraction whereas one with Stot = 2 repels each other, 
which is the origin of magnetization plateaus of the SS 
model |)]]8) . The occurrence of attraction for a singlet- or 
triplet pair is not restricted to the SS system and indeed 
is responsible for the bound states in the 2-leg ladder |l3| ] 
and the S = 1/2 dimerized chain |Ti} |. Interchange of the 
sub (dimer) lattices A— >B ((p x ,p y ) > (p y , —p x )) gives an- 
other hopping matrix corresponding to different bound 
states. Diagonalizing them, we obtain 8 branches (4 for 
each dimer sublattice) for a given value of total spin. Of 
course, on top of them, there are dispersionless bound 
states corresponding to V3 r d and infinitely degenerate 
(up to 5th order) levels corresponding to two isolated 
triplets at co = 2Vb(J, J')- At higher orders of perturba- 
tion, the infinite degeneracy will be lifted and this level 
becomes a narrow continuum. The dispersion curves ob- 
tained this way are shown in Fig. 3-5 for J'/J = 0.5. 
Note that the entire spectrum is invariant under Dad- 
operations {p x ,Py) l-> (Py,Px), {Px,Py) >-> {~Px,Py) 

and their products. 

As is expected, there appear stable branches of bound 
states (bold solid lines and dotted lines) below the two- 
particle threshold (thin broken lines) for S to t = and 1. 
Because of a strong attraction between two triplets on 
adjacent dimer bonds, the energy decrease from the two- 
particle threshold is largest for the singlet bound states 
(Fig. 3) . As a direct consequence of the correlated hop- 
ping starting from (J'/ J) 2 , the dispersion of the bound 
states is relatively large compared with that of the low- 
est single-triplet excitation, which is less than 0.01J for 
J 1 / J = 0.5. To supplement the perturbation theory (PT) 
we also performed exact diagonalization (ED) of the orig- 
inal Hamiltonian for finite clusters with 16 and 24 sites 
under periodic boundary conditions. For small J' / J, e.g. 



0.2, we have confirmed that the results of ED for Stot = 
and 1 are consistent with those of PT. 

We investigated excited states for larger J'/J as well. 
While the bandwidth becomes larger for most branches, 
an almost flat band was found even for relatively large 
J'/ J. For small J'/ J, the energy obtained by ED is 
close to that of the 3rd-neighbor bound pairs (2Vb + V3 r d) 
calculated by PT. Therefore we may identify this flat 
band with the aforementioned 3rd-neighbor bound pairs. 

What is more interesting is that although a quintet pair 
feels repulsion several branches lie slightly below the two- 
particle threshold. This binding beginning at (J'/ J) 3 is 
of purely dynamical origin and is peculiar to the orthog- 
onal dimer systems. Its physical implications to magne- 
tization process are discussed in detail in Ref. Jl7]]. For 
Stot = 2 it is difficult to compare the results of ED with 
those of PT, because there are so many levels around 
the two-particle threshold. However, the minimum en- 
ergy lies below the two-particle threshold as is suggested 
from PT. We examined numerical data for various J'/J 
and found that this seems to be the case at least up to 
J'/J w 0.55, beyond which no decisive conclusion was 
drawn due to finite-size effects. 

Quite recently, Nojiri et al. found in the ESR (elec- 
tron spin resonance) spectra a quintet branch at about 
1400GHz, which is slightly smaller than twice the single 
triplet gap 722GHz We believe that our findings are 
of direct relevance to this observation. 
Selection Rules 

The unique structure of the SS lattice allows only a few 
branches to be observed in INS S. The Fourier trans- 
form S J (p) (j = x,y,z) of local spin operators create 
single-dimer-triplet states 

S^'(p)|G.S.) = / + (p)|j;A(p)) + /_(p)e- i *^+^)|j;B(p)) 

(5) 

over the dimer ground state, where |j;A/B(p)) = 
^a/b j (p) |G.S.) . Since one of the structure factors 

/±(p) = T*sin(Zd(p K ± p y )/2y / 2) vanishes on p x = 
TPy, only an A (B) triplet is excited along the line S: 
Px = Py {Px — —Py)- Although the state S-^p^G.S.) 
contains only a single triplet, the first-order perturba- 
tion (3-point vertices) broadens the wave function and 
S J (p)|G.S.) can have a finite overlap with the two-triplet 
states. Since S^{p x , ±p y )|G.S.) is even under the reflec- 
tion <Td(l,±) (i.e. belongs to Ei-representation), any 
states which are connected to (p x ,±p y )\G.S.) by per- 
turbation should also be even. From this, it follows that 
any bound-state wave functions which contain "6" and 
"c" in a antisymmetric manner (i.e. £2) are orthogo- 
nal to S J (p)|G.S.); that is, INS experiments performed 
along the [110] ([110])) direction (as in Ref. ||) observe 
only bound states shown by solid black lines in Fig. 4. 
Qualitative features agree with what was observed in ex- 
periments [p[. 
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To compare our results with experiments quantita- 
tively, we carried out ED since the value J'/J = 0.635 
recently estimated [|l5| for SrCu2 (1303)2 is too large for 
the results of PT to be used. The results for J = 85K 
are shown in Fig. 6. By close inspection of the wave 
function obtained both by PT and ED, we can identify 
the INS-active branches shown by the lower black lines 
in Fig. 4 in the numerical spectrum for J'/J = 0.635. 
The energies thus obtained are plotted by open circles 
in Fig. 6 giving good agreement with the experimen- 
tal value (filled circles) 0. In particular, the energy of 
the lowest INS-active two-triplet is given by 4.99meV at 
p=0. This agrees with 5.0meV observed in experiments 
H and supports the validity of the value J'/J = 0.635 
for SrCu 2 (B0 3 ) 2 . 

This selection rule does not exclude the possibility of 
observing the remaining branches by other methods. For 
example, optical methods probe excitations at p = 0. 
The representations (E, B 2 ) and (Ai, A 2 , Bi, B 2 , E) are 
active in far infrared spectroscopy (FIR) jl8) and Raman 
scattering respectively. Analysis of the wave func- 
tions shows that, for example, the 1-triplet, the lowest 2- 
triplet (singlet pair), and the second-lowest 2-triplet (3rd- 
neighbor singlet pair) in Fig. 3, 4 belong to E, A 2 © Bi, 
and Ai©B 2 , respectively. This is qualitatively consistent 
with the results of the experiments. Moreover, analytical- 
and numerical results suggest that the lowest bound state 
in the triplet sector is dominated by the 3rd-neighbor 
pair (E-rep), which is not excited by <S J '(p); using the 
above parameters, the energy of it (at p=0) is given by 
36.5cm~ 1 , which does not contradict the FIR- and ESR 
results. Detailed analyses will be published elsewhere. 
Phonons 

Up to now, we have treated dynamics of a purely 
magnetic system. In real materials, coupling between 
phonons and spin degrees of freedom would be impor- 
tant. To see how inclusion of virtual phonons yields ef- 
fective interactions, we treat only the simplest case of 
the Einstein phonon where each dimer (J) bond rotates 
independently around its equilibrium position with a fre- 
quency loq. In contrast to the purely magnetic case, pair 
creation/annihilation and one-particle hopping which ac- 
company emission of a phonon are allowed. 

Here we only consider the most important case of a 
nearest-neighbor pair. A little algebra shows that the 
second-order processes generate an attractive interaction 
V^f 1 cx — (2J + ujq)/(J + ^0) even for a quintet pair; 
a singlet pair feels attraction while a triplet repulsion. 
Therefore, we may expect that the lowering in energy for 
quintet bound states is enhanced by phonons. 

In SrCu 2 (B03) 2 the interlayer coupling J± is not neg- 
ligible either. However, J± changes nothing as far as 
inira-layer excitations are considered Jl5| ] . It might make 
a new (INS-inactive) interlayer bound states slightly be- 
low the two-particle threshold by a very weak interaction 



Jl(1 - l/4(J'/J) 2 )T r T 2 . These bound states have a 
small bandwidth of the order ((J'/J) 6 ). 

In conclusion, we demonstrated that unique dynamical 
property of the SS model allows an interesting pairing of 
triplets. Bound two triplets can move on a lattice much 
easier than a single isolated triplet can and this fact ex- 
plains the recent INS experiments. Note that the differ- 
ence in dispersion is not quantitative and originates from 
that in hopping mechanisms. Moreover, it enables the 
two triplets to form quintet 'bound' states despite the re- 
pulsion between them. To complement the perturbative 
consideration, we also carried out exact diagonalization 
and obtained for J'/J = 0.635 the results consistent with 
those of experiments. 

After this work was comp leted, we became a ware of 
a preprint by Knetter et al ( cond-mat 0005322 ). They 
investigated gaps at p=0 for singlet- and triplet bound 
states using a similar but different method. 
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CAPTIONS 

Fig. 1 

Configurations of two orthogonal dimers in a unit cell. If 
the direction of the vertical {a) bond is fixed (say, point- 
ing downward), two configurations (left and right) are 
allowed for the horizontal (/?) bond. 
Fig. 2 

"Hopping" processes of two triplets. A(B) dimers are 
shown by bold black (gray) lines. A filled circle denotes 
a representative point of a unit cell r. Symbols for the 



matrix elements are the same as those in Ref. 17 . 
Fig. 3 

Singlet (5 to t = 0) dispersion in the [110]- and [100]- 
direction for J' /J = 0.5 obtained by the perturbation. 
Dashed- and dotted lines denote B-branches and 3rd- 
neighbor bound states, respectively. Two-particle thresh- 
old lies at 1.356 J (thin broken line). 
Fig. 4 

Same for triplet (<Stot = 1) sector. Only branches shown 
by solid black lines are observable in neutron-scattering 
experiments. 
Fig. 5 

Same for quintet (5tot = 2) sector. Note that several 
branches lie below 2-particle threshold (thin broken line) . 
Because of higher-order processes, states above the 2- 
particle threshold (thin dashed line) are unstable. 
Fig. 6 

The branches with S t ot = 1- The results of the ED for 
J' /J — 0.635 and J = 85 K are shown by open sym- 
bols. Among them the branches observable in neutron- 
scattering experiments are shown by the open circles 
(lower: 1-triplet, upper: 2-triplet). For comparison, 
the experimental results are shown by the closed circles. 
Small splittings at p=0 are artifacts coming from the 
shape of a cluster (which is ^-invariant) used in ED. 
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Fig. 1 Totsuka, Miyahara, and Ueda 




Fig. 2 Totsuka, Miyahara, and Ueda 
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Fig. 3 Totsuka, Miyahara, and Ueda 
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Fig. 4 Totsuka, Miyahara, and Ueda 
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Fig. 5 Totsuka, Miyahara, and Ueda 




Fig. 6 Totsuka, Miyahara, and Ueda 



